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A sophisticated analysis of EGRET data has found evidence for gamma-ray emission from the 
galactic halo. I entertain the possibility that part of the EGRET signal is due to WIMP annihilations 
in the halo. I show that a viable candidate with the required properties exists in a model with an 
extended Higgs sector. The candidate has a mass of 2-4 GeV, a relic density Q. ~ 0.1 (for a Hubble 
constant of 60 km/s/Mpc), and a scattering cross section off nucleons in the range 10~^-10~^ pb. 
The model satisfies present observational and experimental constraints, and makes strict predictions 
on the gamma-ray spectrum of the halo emission. 



Dixon et al. [Q, in a sophisticated analysis of the 
EGRET gamma-ray data, have found evidence for 
gamma-ray emission from the galactic halo. Filtering 
the data with a wavelet expansion, they have produced a 
map of the intensity distribution of the residual gamma- 
ray emission after subtraction of an isotropic extragalac- 
tic component and of expected contributions from cosmic 
ray interactions with the interstellar gas and from inverse 
Compton of ambient photons by cosmic ray electrons. 
Besides a few "point" sources, they find an excess in the 
central region extending somewhat North of the galactic 
plane, and a weaker emission from regions in the galactic 
halo. They mention an astrophysical interpretation for 
this halo emission: inverse Compton by cosmic ray elec- 
trons distributed on larger scales than those commonly 
discussed and with anomalously hard energy spectrum. 

I find it intriguing that the angular distribution of 
the halo emission resembles that expected from pair an- 
nihilation of dark matter WIMPs in the galactic halo, 
and moreover, that the gamma-ray intensity is similar 
to that expected from annihilations of a thermal relic 
with present mass density 0.1-0.2 of the critical density. 
Namely, the emission at 6 > 20° is approximately con- 
stant at a given angular distance from the galactic center 
- except for a region around {h,l) — (60°, 45°), corre- 
lated to the position of the Moon, and a region around 
(&, I) = (190°, —30°), where there is a local cloud (6 and 
I are the galactic latitude and longitude, respectively). 
Dixon et al. [|l| argue against the possibility of WIMP 
annihilations on the base that direct annihilation of neu- 
tralinos into photons would give too low a gamma-ray 
signal. However most of the photons from WIMP an- 
nihilations are usually not produced directly but come 
from the decay of neutral pions generated in the particle 
cascades following annihilation. 

In this letter I show that WIMP annihilation can ac- 
count for both the intensity and the spatial distribution 
of the halo emission. A model independent analysis, al- 
though instructive (see Ref. |0]), would necessarily be 



vague. So to be concrete, I introduce a particle physics 
model with a suitable dark matter candidate that satis- 
fies present observational and experimental constraints. 

Consider a few GeV Majorana fermion in a model with 
an extended Higgs sector. The Higgs sector contains two 
Higgs doublets Hi and H2 and a Higgs singlet N . Let 
the Higgs potential be 



l^tiggs = Al (hIHi - vf] + A2 (hIH2 - V. 
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Taking all the A's real and positive guarantees that the 
absolute minimum of the potential is at (Hi) = wi, 
(H2) — V2, and {N) = vn, with zero vacuum expectation 
values for all the other fields. 

There are 6 physical Higgs bosons: one charged Higgs 
boson H^, two neutral "pseudoscalars" Pi and P2, and 
three neutral "scalars" Si, S2, and 5*3 (I label the Ps 
and Ss in order of increasing mass). The charged Higgs 
boson is ~ iJ~*sin/3 + H2 cosf], and its mass is 
= Xq^^v, where tan/3 = V2/V1 and v = y^vl + f|. 

The two pseudoscalars are linear combinations Pk = 
U[i A + [7^2 iV/ of A = V2{sin (3 Im Hf + cos /? Im ) and 
Nj — \/2ImiV. The unitary matrix diagonalizes 
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The three scalars are linear combinations S 
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U^2H2R + Uf^NR of HiR = \/2 Re iJ?, H2B. = \/2 Re H^, 
and Nr — \/2Re A^. The unitary matrix diagonahzes 



where A, 



4AiVi + A45UI 
A45U1W2 
2XzV2Vn 

i = Xi + Xj . 
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The Higgs bosons couple to the ordinary leptons and 
quarks and to the new Majorana fermion x through the 
Yukawa terms 

-Cvukawa = hdQHidR + huQH2UR + hiLHiBR + h^xxN. 

(4) 

Here Q — {uL,dL) denotes the SU(2) quark doublets, 
L = the SU(2) lepton doublets, and dji,UB,eii 

the SU(2) singlets. Generation indices are suppressed. 

When the Higgs fields get vacuum expectation val- 
ues, the up-type quarks, the down-type quarks and the 
charged leptons acquire masses to„ = huV2, rrid = hdVi 
and mi — hivi, and the new Majorana fermion acquires 
mass — hvjy. The interaction terms of the Higgs 
fields with the quarks, the leptons and the x read: 
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Now I find the interesting region in parameter space 
in which the calculated gamma-ray emission matches the 
Dixon et al. maps while satisfying present observational 
and experimental constraints. 

The gamma-ray intensity from xx annihilations in the 
galactic halo is given by 



4>^{b,l,E) = n^{E) 



47rm^ 



p\dl, 



(6) 



where (j)^{b,l,E) is in photons/ (cm^ s sr), n^{E) is the 
number of photons per annihilation with energy above 
E, av is the xX annihilation cross section times relative 
velocity at v = 0, py- is the x mass density in the halo, 
and the integral is along the line of sight. For a canonical 
halo, p{r) — p\oc{r1 + R^) / {r1 the integral depends 

only on the angular distance "0 from the galactic center, 
and is practically independent of the halo core radius Tc 
a.t ip — 56°, where it is / p^dl ~ 2.8p\ocR- Notice that 
for this halo model, the intensities a,t ip = 40° and at 
"0 = 60° are approximately in the ratio 2:1 as on the 
Dixon et al. maps. 

I consider interesting those parameter values for which 
the calculated gamma-ray fiux at = 56° is between 
7 and 8 times 10^^ photons/cm^ /s/sr for a local dark 
matter density pioc in the range 0.3-0.5 GeV/cm^ and 
a distance from the galactic center R = 8 kpc. Fig. 1 
shows the interesting region in the av-m^ plane (a point 
for each choice of parameter values). The figure includes 
the bounds discussed in the following. 




[GeV] 



FIG. 1. Region in the av 



plane where the x particle 



can explain the gamma-ray halo signal and satisfy all experi- 
mental constraints considered. 

The most important constraint on the interesting re- 
gion comes from the measured fiux of cosmic ray antipro- 
tons. If the gamma-rays are produced in jets originated 
by quarks, there is an associated production of antipro- 
tons. The ratio of antiproton and gamma-ray fiuxes is 
independent of the WIMP annihilation cross section and 
of the local mass density, and the relative number of an- 
tiprotons and photons per annihilation is fixed by the 
physics of jets to a number of order 1. Since the antipro- 
tons are confined in the galaxy for ~ 10^ yr while the 
gamma-rays from the halo take no more than ~ 10^ yr 
to reach us, the antiproton flux is ~ 10"^ times larger than 
the gamma-ray flux. More precisely, the antiproton flux 
at a p kinetic energy T is 



, dNp av 2 , 



(7) 



For the antiproton spectrum per annihilation dNp/dT 
I use the Lund Monte-Carlo, and I evaluate the con- 
tainment time fcont in a diffusion model jsj: tcont — 
(1 + p/GeV)-°-65 X 10^5 s. The factor p = [T{T + 
2TOp)]/[(T + A){T + A + 2mp)] takes into account so- 
lar modulation. Given the uncertainties in the antipro- 
ton propagation in the galaxy and in the effect of the 
solar modulation, I accept models with 0j3(2OOMeV) < 
3 X 10~s p/cmVs/sr/GeV g (taking A = 600 MeV). 
This bound effectively limits the viable decay channels 

to XX T+T-. 

The X relic density can be obtained with a standard 
procedure S, 



1.023 X 10"27(jjj^3/g 
Jq^ {<jv)gl^^dx 



(8) 



where Xf is the solution of Xj:^ + ^ln{g^/xf) = 80.4 -|- 
ln{m^(crw)) with m-)^ in GeV and (av) in cm^/s. Here 
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{av) is the thermally averaged annihilation cross section 
at temperature xm. Since freeze-out occurs below the 
QCD phase transition, I take 51* = 81 for the relativistic 
degrees of freedom at freeze-out. The result of the relic 
density calculation is shown in Fig. 2. The relic abun- 
dance in the interesting region turns out to be in the 
cosmologically interesting range: Vl^ ^ Q.l for a Hubble 
constant of 60 km/s/Mpc. 
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FIG. 2. The x relic density for x particles that can ex- 
plain the gamma-ray halo signal and satisfy all experimental 
constraints considered. 

In the interesting region, two Higgs bosons, 5*1 and 
Pi, are light, with masses between 2 and 16 GeV. They 
are mostly singlets and their production in colliders is 
very suppressed. First, I consider the LEP bounds 
on the search of Higgs bosons in two-doublet models. 
Since there is an additional Higgs singlet, sin^(/? — a) 
in the bound from e+e^ Zh must be replaced with 
\Uii cos f3+Ui2 sin and cos^(/?— a) in the bound from 
e+e~ hA must be replaced with \U12U12 - J7nC/fip. 
In the allowed region, the reinterpreted s\T?{f3 — a) and 
cos^(/3 — a) are smaller than 5 x 10~^, and are not ex- 
cluded by accelerator searches Secondly, I consider 
Higgs bremsstrahlung from final state leptons and quarks 
in Z decays. I find 

I{Z ^ f-JX) ^Gf (mx\ 2,2 



^ ,//) 



mz 



where Ap.f — Ul{ col j3, As^f — [//^/cos/? for up-type 
quarks, Ap.j — U^ta.n(3, As^f — C/j|/sin/3 for down- 
type quarks and leptons, c/ = 3 for quarks and c/ = 1 
for leptons. The function g{y) comes from the phase- 
space integration, and is g{y) = 1 at mx = 12 GeV. The 
most stringent constraint from the LEP results on two 
leptons -I- two jets production is T{Z ffX)/T{Z 
fj) < 1.5 X 10--* at mx = 12 GeV. The highest values 
I find in the interesting region of my model are of order 



10~^, and so this constraint is easily satisfied. 

In addition to accelerator bounds, the x particle must 
satisfy present bounds from direct and indirect dark mat- 
ter searches. The cross section for elastic x-proton scat- 
tering, which is spin-independent, is 



o 4 4 



{nip 




where kd = {mddd + russs + nibbb) — 0.21 and fc„ = 
{myllu + rricCC + m-ttt) = 0.15. Fig. 3 shows a compar- 
ison of this cross section with the presently most strin- 
gent upper bound on few GeV WIMPs 0. The bound 
touches the interesting region, but to explore it all major 
improvements in sensitivity are needed, especially at low 
thresholds f§. 




FIG. 3. The x~proton spin-independent scattering cross 
section versus the x mass. The solid line is the present bound 
from ref. [7]; the dots show the interesting region. 

The X fermions can also accumulate in the Sun and in 
the Earth cores, annihilate therein and produce GeV neu- 
trinos. Accumulation in the Earth is not efficient because 
the X mass is below the evaporation mass from the Earth 
(~ 12 GeV), and no neutrino signal is expected from 
there. The evaporation mass in the Sun 3 GeV) is in 
the middle of the interesting region. Following Ref. , I 
have calculated the neutrino flux from the Sun including 
evaporation, annihilation and capture, and then the rate 
of contained events and the flux of through-going muons 
in neutrino telescopes. I have imposed the experimen- 
tal bounds in Ref. |jl^. Fig. 4 shows that, although the 
present limits on through-going muons constrain the al- 
lowed region, it will be difficult to probe models with 
below the evaporation mass of ~ 3 GeV because of the 
exponential suppression of the signal. 
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FIG. 4. The flux of through-going muons from the Sun 
induced by neutrinos from x annihilations versus the x mass. 
The solid lines are the Baksan and MACRO limits [10]; the 
dots show the interesting region. 

The clearest test of WIMP annihilations in the halo is 
the shape of the gamma-ray spectrum. It must show 
a sharp cutoff and a gamma-ray line at the WIMP 
mass. The continuum spectrum is determined from 
jet physics. The gamma-ray line intensity is propor- 
tional to the annihilation cross section into two photons, 
which in the specific model presented here results in the 
range 10~^^-10~'^'' cw?/s (using formulas in Ref. pH). 
The poor energy resolution of present gamma-ray de- 
tectors unfortunately smears the line out: even a tan- 
talizing energy resolution of 0.1% merely gives 10~^^ 
photons/cm^ /s/sr/MeV in this model. Fig. 5 illustrates 
the gamma-ray spectrum expected from xx annihilations 
for X particles in the interesting region. For comparison, 
the figure indicates the measured extragalactic flux . 
Since the antiproton constraint limits the x mass to less 
than 4.3 GeV, a clear detection of halo gamma-rays more 
energetic than this would be an indication against the 
present WIMP model. 

In conclusion, it is possible to explain the Dixon et 
al. gamma-ray halo as due to WIMP annihilations in the 
galactic halo. As illustrated in an explicit model, not only 
the intensity and spatial pattern of the halo emission can 
be matched but also the relic density of the candidate 
WIMP can be in the cosmologically interesting domain. 
The model is quite predictive and can be tested with a 
variety of techniques, above all by accurately measuring 
the energy spectrum of the halo gamma-ray emission. 

I thank D. D. Dixon for discussions. 
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FIG. 5. The predicted gamma-ray spectrum from WIMP 
annihilations in the halo compared with the measured extra- 
galactic flux [12]. Solid lines correspond to x masses of 2.3, 3, 
and 4.3 GeV, and they delineate the range of expected signals. 
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